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Summary 


This report presents the full Navier-Stokes time- dependent , compressible, 
turbulent mean-flow equations in mass-averaged variables for plane or axi symmetric 
flow configurations. The equations are derived in a body-oriented orthogonal 
curvilinear coordinate system. Turbulence is modelled by a system of two equa- 
tions for the mass-averaged turbulent kinetic energy and dissipation rate which 
determine the turbulent (eddy) diffusivity, recently proposed by Wilcox et al. 
These equations are rederived for the coordinate system and flow configuration 
considered, and some new features of these equations are discussed. A system of 
second-order boundary layer equations is then derived which includes in a con- 
sistent way the effects of longitudinal curvature and the corresponding normal 
pressure gradient. In this system the normal momentum equation is retained. The 
Wilcox and Chambers approach is used in considering effects of streamline curva-- 
ture on turbulence phenomena in turbulent boundary layer type flows. Their two- 
equation turbulence model with curvature terms are rederived for the cases 
considered in the present report. The derived system of seven second-order 
boundary layer equations serve as a basis for an analytical-numerical investiga- 
tion of a variety of boundary layer (parabolic) type problems where streamline 
curvature is of the order of the characteristic length In the longitudinal 
direction. 
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speed of sound 

specific heat at constant pressure 

mass- averaged turbulent kinetic energy 

total (stagnation) enthalpy 

specific static enthalpy of the mixture, h = 

specific static enthalpy of species i 3 

diffusional flux factor 

heat conductivity 

characteristic length of flow in x-direction 

Lewis number 

dissipation length scale 

molecular weight 

Prandtl number 

pressure 

heat conduction vector 

universal gas constant; radius of curvature of the body surface 

Reynolds number 

distance to x-axis 
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x- component of velocity 
y- component of velocity 
total velocity vector 

net mass rate of production of chemical species 
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longitudinal coordinate along surface of the body 
coordinate normal to the body 

species mass fraction; angle of body with x-axis; turbulence 
closure coefficients 
turbulence closure coefficients 

specific heat ratio 

characteristic length of flow in y- direction 

mass-averaged turbulent dissipation rate 

curvature of the body 

coefficient of viscosity 

kinematic viscosity 

eddy diffusivity 

turbulence closure coefficients 

density 


„ turbulence closure coefficients 

c r* 

'W viscous stress tensor 

cp cylindrical or azimuthal coordinate 

gj vorticity or pseudovorticity 

Q dissipation rate quantity, poo 

^ divergence operator 


Subscripts 

a a = 0 plane flow, a = 1 axi symmetric flow 

t denotes turbulent quantity 

i denotes species 
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Superscripts 


~ mss- averaged value 

* nondimensional quantity 

— time- averaged value 

f time- averaged fluctuations 

" mass-averaged fluctuations 



1 . 


INTRODUCTION 


Recently, considerable research efforts have been devoted to the 
understanding of the aerodynamics of combustion and the complex turbulent 
mixing phenomena occurring in combustion devices- Experimental investigations 
and theoretical modelling of such phenomena combined with sophisticated computa- 
tional fluid- dyn a mic methods are the essential tools facilitating the economical 
design and operation of these combustors. 

Usually all practical flames have the form of turbulent jets issuing 
from round orifices, the fuel gas being introduced through a central jet and 
the oxidizer (air) through an annulus surrounding it, so that straight concentric 
jets are formed. However, in order to enhance the fuel-air mixing process, the 
primary (fuel) jet or the secondary jet (air) or both are given a certain degree 
of curvature (e.g., swirl). Also, in certain types of combustors involving 
coiipressible jets, in order to achieve optimum combustion characteristics, the 
finite-rate chemical processes should be controlled by controlling the pressure 
and temperature fields in the combustor. This can be achieved by imparting to 
the flow various degrees of curvature. 

Among the theoretical models describing curvature effects on turbulence, 
second-order closure two-equation turbulence models, which utilize two parameters 
to characterize the turbulence and to determine the eddy diffusivity, with each 
parameter satisfying a nonlinear transport equation, involve less empiricism 
than the mixing-length theories or the one equation models of turbulence. They 
require no advance knowledge of the flow under consideration. Yet they are 
simple enough to use for general engineering applications. Such models have 
proliferated. Recently, Wilcox and Chambers (Ref. 2) have proposed an 
interesting extension of their model to compressible flows with streamline 
curvature. The success achieved by this method in predicting the behaviour 
of a number of turbulent boundary- layer flows (Ref. 2) warrants its further 
use and testing in other compressible applications. 

The present report is concerned with the analytical investigation of 
curvature and compressibility effects on the turbulent mixing of a compressible 
(fuel) jet issuing into a compressible air flow (Rig. l) . Flow variables of 
primary interest are the fuel/ air ratio and temperature distribution. The 
full Navier-Stokes equations for global continuity, species continuity, 
momentum and energy conservation are derived for the mean mass -averaged varia- 
bles of a compressible, turbulent, axisymmetric or plane flow in curvilinear 
orthogonal body-oriented system of coordinates (Ref. l) . The apparent (Reynolds) 
stresses were modelled according to the eddy viscosity concept, the turbulent 
viscosity being a function of two parameters, the turbulent mixing energy and 
the turbulent- dissipation rate determined by two nonlinear partial differential 
equations. The Navier-Stokes form of these equations was obtained in the 
coordinate system considered, with due account of compressibility and curvature 
effects. A compressibility term was found to exist in the equation for 
turbulent dissipation omitted by Wilcox and Chambers in their general equations 
in Cartesian coordinates. Starting from this system of seven equations, 
a complete set of second-order boundary-layer (parabolic) type equations was 
derived (terms of order unity and of order 5t/R are retained, 5^ "being the 
characteristic mixing zone thickness and R the radius of curvature, Which is 
assumed to be of the same order of magnitude as the characteristic longitudinal 
length of the problem) . 
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2. MASS -AVERAGED EQUATIONS OP MOTION 

The general equations describing flows of a reacting mixture of 
perfect gases in vector form are given by: 


Global Continuity: 


- P r-° 


Species Continuity: 


&(p « ± ) 


at 


+ f • (pOjV* + 3t) = w ± 


Momentum: 


^ • ptftf = ~ ^ ^ 


Energy: 


EpHzE + y 

at 


pV*H + q +^V - & ■ V) 


= 0 


Equation of State: 


— i OC . 

P = HTP > 


1 — i 


( 1 ) 

( 2 ) 

(3) 

( 4 ) 

(5) 


where ^ is the divergence operator; p, T, p are the density, temperature and 
pressure of the mixture, respectively; v the velocity vector; aq the mass 
fraction of the i^h species; jj_ the diffusional mass-flux vector, jq = p(VjL“V); 
Wq the net mass rate of production of species i per unit volume by chemical 


reaction (grams of i per cm3 per sec); ^ the viscous stress tensor; H = ? ah 

+ f • 2 the total (stagnation) enthalpy; hq the specific enthalpy of the 

i'tih species; cf the heat- conduction vector; Mq the molecular weight of the i^* 1 
species; R the universal gas constant. Constitutive relationships for the 
viscous stress tensor ^ , for the diffusional mass flux vector j , and for the 
heat conduction vector q* must be added to these equations in order to obtain 
a closed system. Thus, for Newtonian fluids, i.e. fluids such that the 
viscous stress tensor is a linear function of the rate of strain, the viscous 
stress is given by: 



def t- 


2 

3 


div 


) 


( 6 ) 


2 



with def f = grad V* + (grad the superscript T denoting the transpose of 

a tensor 3 p. being the coefficient of viscosity and X the unit tensor. The 
diffusion of each species is assumed to depend only on the gradient of the 
particular species mass fraction. The assumption requires that the binary 
diffusion coefficients Dy for each species are equal and the simple Fick*s 
law is applicable. Thus 

3*1 = - PDy grad ^ H ff Srad CC ± (7) 

where Le and Pr are the Lewis and Prandtl numbers, respectively. |^irthermore , 
the fluid is assumed to obey Fourier* s law of heat conduction for q: 

c fi 

q* = - k grad T = - p^- grad T (8) 

where Cp = § <X^_Op^ 9 Cp^ being the specific heat of species i. 

For turbulent flows considered in the present investigation, it is 
assumed that the functional form describing the transport processes discussed 
above remain unchanged and the laminar coefficients are replaced by their 
turbulent counterparts, i.e. p. etc. The coefficients ... have yet 

to be determined. 

The body-oriented orthogonal curvilinear coordinate system used is 
shown in Fig. 1, where x is the distance along the surface of the body (measured 
from a certain point 0) and y the distance normal to this surface. The corre- 
sponding components of the velocity vector V* are u and v. R(x) is the radius 
of curvature of the body, reckoned positive for a convex body, and a(x) its 
angle with and r-b(x) its distance from the axis. The length element is given 
hy 


jn 2 , 2 , 2 ^ , 2 J 2 L . 2, 2 

di = h^ dx + h 2 dy + h^ dcp (9) 

where the metric coefficients are 

h^ = 1 + Ky, h 2 = 1, h^ = r a , r = r^ + yco£#, k(x) = l/R(x) (10) 

a = 0 for plane and a = 1 for axi symmetric flows, and cp is the lateral (z.) 
coordinate in plane flow and the azimuthal angle in axi symmetric flow. Then 

from the usual relations for vector operators in orthogonal curvilinear coor- 

dinates the continuity equation, Eq. (l), is found to be 

it + > ( ^ + ; y) - { s ( papu > + b } = 0 (u) 

Using Eq. (7) 5 the species continuity equation, Eq. (2), becomes: 


3 
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r a (l+Ky) 


5a. -1 
Sc dy 


+ W, (12) 


where Sc = Pr/Le is the Schmidt number. 

The Navier-Stokes momentum equations are (Ref. 3): 


x-momentum: 


^ + ~ a (l+Ky) ife [Ap^Vt^)] + [r a (l+Ky)(puw-T xy )] j 


+ ifey (p^-V ‘ r (v~ T w> silCl = 0 


(13) 


y- momentum: 

d (pv) + L 


"ST 


r a (l+Ky) 


|e + b [^(^^(P^V^yy)] 


( pu +p— T ) " *- (p-T )cOS0£ = 0 

1+Ky K * xx y r x qxp y 


(i4) 


The viscous stresses t-^-, . . . > t ^ are given by Eq. (6), and in the coordinate 
system considered are: 


T = LI 
XX ^ 


2 / du 


1+Ky V 5x 


+ KV 


I div V* 


T = [l 
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^ Sv 2 rr> 

2 S - 3 T 


T = |I 
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t = [1 
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div 
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With Eqs. (7-8) and relations: 
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5E “X h i (T) + X a i °Pi i§ + |e (?" 


If = I h i (T) 3T + Z a i °Pi § + |f (r n ' 
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where Cp. = dh^/dl, V 2 = (u 2 +v 2 )/2 J and assuming Le = 1, the energy equation, 
Eq. (4) , 1 becomes: 


^ + 7^{^ + 5 [ra(1+Ky)wH1 


(l+Ky) 


r a (l+Ky) 


r a 

r a [x 5 f V 2 \ 

. a 

15E 

l+Ky Pr Sc \ 2 J 
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+ 5E t r ^ T xx u+T xy V)] + Ip [r a (l+Kry)(T xy u + T yy v ) ]} 


■f): 

(16) 


An alternate form of the energy equation may be obtained if we make 
use of the equation for the kinetic energy of the flow. Multiplying Eq. (13) 
by u and Eq. (l4) by v and adding, we obtain, after some calculation: 



— - v 

1+^y dx 
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U+T xy v)] + 


a 

3y 
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XX 
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(17) 


Noting that H-(V^/2)=h and using Eq. (17) , we derive from Eq. (l6) the energy 
equation in the following second form: 


^tr 1 + 7 a ( - ~ — { s (r&puh) + If [* a (i+Ky)pi*]} = 


1 f f r a \x dh \ d 
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u 5p 
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(18) 

Contd. • . 
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In order to obtain the governing conservation equations for turbulent 
flows , it is convenient to replace the instantaneous quantities in the above 
equations by their mean and their fluctuating quantities. In this treatment, 
the mass- weigh ted- averaging procedure is used. With this procedure the 
resulting- equations for the average turbulent flow quantities will have a 
form very similar to the equations for laminar flow. 

If f is any flow quantity, the conventional time average of this 
quantity is denoted by f . Then, the mass average of any quantity except 
density, pressure and viscous stresses (fluctuations of transport coefficients 
are neglected) is given by 


f 


The quantity f may then be written as 

f = f + f" 

where f" is the departure or fluctuation from the mass-averated value f . It 
shou ld be noted that f M ^ 0; f" = -p'f T /p. Also, it is easy to show that 
pf IT = 0. Representing the density, pressure and viscous stresses by the sum 
of their time average and its departure or fluctuation from the time average, 
namely 


= f 

p 


P = P + P 1 » P = P + P' s T ij = T ij + T ij (19) 

and the velocities, the total enthalpy, static enthalpy, temperature and mass 
fractions by their mass average and its departure from the mass average, 
namely 


u = u + u", v = v + v", H = H + H", h = h + h" 

~ ~ ~ pa. p. \ 

T = T + T ,: , a. = a . + aV ( a. = — - = — ) 

5 i x i V i - / 

v p p / 

substituting these quantities into Eqs. (ll)-(l8), and applying the Reynolds 
rules of averaging (Ref. 4) yields the following conservation equations for 
plane or axisymmetric turbulent flows: 


6 



Global Continuity: 


I + 7^ { 5E <'*® + 1 I } - o (a) 


Species Continuity: 


+ _^ 

dt r a (l+Ky) 
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2 r a 
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+ t^ 


r a ( tL _L_ . paTEF 
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(22 > 


x-momentum: 


a , a r„a /n t! k'pot = 1_ dp + 


^ + 7^ 1 S + 5 tr-(l*x)p5r^ * S| - - 1 


r a (l+ny) 
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[r (t vv - pu"u")] + ^ [r (l+Ky)(T vv - pu"v")] f + 


3x L r ^ T xx 
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(23) 


y-momentum: 
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(24) 


Energy: 
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(25) 

Contd. 
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h { r& f h H + (; xx + T ix^ + u ") + ( V + T iy )( ^ + v,,) ] } + 


+ If { ra d + ^) [ fe § + ( V + T iy )(S + u " } + ( V + T ^ )( ^ + v,,) 
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+ U pu V + V pv V 


/ ~2 t 
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where 


pu ,T u M + py ,, v ,T 


is the mass averaged specific turbulent kinetic energy. Hence Eq. (25) can 
be rewritten as 


It [r a (l + K y)(p H - p)] + 4- r a p u ( h + — - * — ■ - + e Xl + 
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(27) 

Contd. 
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(27) 


+ h{ ra £* ( ^oc - + ^ v - 

+ ly { r& ( 1 + ' pu " v ") + ^ V " pv^]} 

In Eqs. (22), (23), (24) and (27) the Reynolds (apparent) stresses 
and the turbulent mass and energy transports are calculated by applying an 
eddy transport concept which is postulated as: (Ref. 6) 
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where vt is the turbulent (eddy) diffusivity, and Pr-fc and Set the turbulent 
Prandtl and Schmidt numbers, respectively. It should be noted that unlike 
I a, the molecular viscosity, vt is not a property of the fluid. Its value 
varies from point to point in the flow and is determined by the structure 
of turbulence at the point considered. Substituting Eq. (28) into the 
conservation equations, Eqs. (22)-(24) and (27) we get 
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x-momentum: 
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Energy: 
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— — — . r>* / “^2 ~2\ / 

where P = p + 2pe/3 ? H = h + (u + v )/2 and the turbulent kinetic energy 
has been neglected compared to the mean kinetic energy for a high-speed 
flow. 


Using the equation for the mean kinetic energy of the turbulent 

flow, pV 2 /2 = p(u 2 + v 2 )/2, obtained by multiplying Eg. ( 30 ) by u and Eg. 
( 31 ) by v and adding 
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ef f , , - \ ( ^ usincc + vcosa 2 . ^ 2 - 

= (|i + pv t ) ( 2a - div V J - ^ pe 


eff t \ ( 1 Sv'.Bu k ~ 
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Substituting the above into Eq. (32) we obtain the following alternate form 
of the energy equation: 
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Equations (21), (29), (30), (3l) a (32) or (34) and (5) form the 
complete set of Navier-Stokes equations describing plane or axi symmetric 
laminar and turbulent reacting flows containing i distinct species. The 
turbulent diffusivity v-fc is yet to be determined in terms of known or 
calculable quantities. Following Wilcox and Traci (Ref. 5) it is assumed 
that the eddy diffusivity is expressed by the equation 



(35) 


2 

where e is the mass-averaged dissipation rate and is defined as e = pft H ft"/p 
where ft = ft + ft" = poo, 00 being the vorticity or pseudo-vorticity . Both e and 
e are determined from the following two partial differential equations (for 
a detailed derivation of the equations for the particular coordinate system 
considered here see Appendix A) . 
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M££l + _ 

£yfc „a 


r*(krt { S (r "° S e) + b [ra(l + W) ' P " e] } * 


= (of* p S - £* e)e - e p div V + 


b r r ^ + °* py &e 

r a (l+ K y) 1 S 


— {' 

f 1 +ky) L 


1 + Ky Sc 


+ b 


\l + k3T)(m. + cr* pv t ) || J 


(36) 


and 


S(p S) , 1 f S / a - ^ 2^ S r a /-. n- ~ 2 n l 

+ a" — {S (r »“'l + 5 [r(l*.j)pve)}- 


r (l+Ky) 


o; p s e 2 - -I p + 2(7 


i a* Y . f ae 

l+Ky Sc J 


2 -i 


e 3 + 


r a (l+Ky) 


a 

SE 


r a (tx + crpv t ) ^2 ^ 


l + Ky &c 


cSy 


r a (l + Ky)(p. + apv t ) j-- 


- 2 ~ 
p e div V 


» pe div V / dp dp _ dp dp 
S 3 ~2 1+ny V Sc dy ” ^y Sc 

a x 


(37) 


The term in box 1 in Eq. (37) is absent in the corresponding equation given in 
Ref. 5 by Wilcox and Traci for general three-dimensional Cartesian coordinates. 
The term in box 2 is neglected in Ref. 5; here it is kept and modelled (see 
Appendix A for details). The importance of these terms* as well as the values 
of the closure constants I2 and £3 can be assessed by comparison with appropriate 
and reliable experimental results. The turbulent Prandtl and Schmidt numbers 
Pr£ and Set in Eqs . (29) and (32) or (3*0 and the closure coefficients a* a*, 

P* p** cr, cr** g* appearing in Eqs. (35) and (36) are (see Ref. 5 for a detailed 
discussion) : 


P = 


_3 
20 ? 


P* 


9 

100 5 





(38) 

Conti. . . 
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1 " II exp (~ Re t/ 2 ) 


■ °* = IS 


1 “ II e2iI) (" 2 Re t) 


i = pe^/e, Re t = pe 1 / 2 i , Pr t = Sc t = | (Le = l) 


( 38 ) 


j0 being defined as the length scale of turbulence 


S = \l2(e i e e + 2e e Y 
v xx xx yy yy xy yx ; 


(39) 


l 5u ~ ^ 5v 

e xx = 1+^ V S + KV / e yy = 


= 1 / _1 dv <Ju _ ku 

e xy 2 V 1+ K y ^x ciy ~ 1+ K y 


(kO) 


To complete the formulation of the set of equations an appropriate 
set of boundary conditions must be specified. Generally, for nonturbulent 
regions e = e = v-fc = 0, while for solid boundaries the usual no-slip boundary 
condition on u and v applies. Moreover, either the temperature or heat flux 
must be specified. The turbulent energy and length scale must satisfy the 
condition e = Jl = 0 (see Ref. 5) • 


3. E OUKDARY - LAYER TYPE EQUATIONS OF MOTION 


To compare the relative magnitude of the terms in the equations 
of motion for the mean turbulent flow given in the previous section in 
situations where the characteristic length of the flow in the x direction 
is L and in the y direction is &fc « L, we introduce the following new dimen- 
sionless variables that are of order unity in the flow field considered: 


R 


U t 
e 


X* = T > K* = r* = y* = t* = 


u* = v* = p* = p* = 


_ p_ 


* 


* w 


e* = 


2 5 


e6. 


e* = 


H* = 


P U 

e 


U 


PeV 


n*. - * - P^e* l/2 


(4l) 


Ik 



where U e is a characteristic velocity in x-direction, p e a characteristic 
density and and Y’t can be an y velocity characteristic of the turbulent 
field. The characteristic velocity V e is yet to be determined. 


In nondimensional form the continuity, Eq (21) , is: 

1 


dp* . 

5P 


r* 


1 + 


L H *) 


|p: ( r* a p* u*) + 


L V e d 

6 t U e^ 


r*" 


1 + — ^ ) p* v* 
L R* J ^ 


(42) 


If all the terms of this equation are to be of the same order of magnitude 
then 


L V e 
6 t U e 


and 


U 5, 
V =-£-t 
e L 


(43) 


Thus the continuity equation, Eq. (42) , becomes 


dp* 

dt* 


r*° 


i 

L R* 


{ is? p* s *) + 


r* 


1+ r§)"p*^ 


}■ 


(44) 


The species continuity equation, Eq. (29), in nondimensional form, is 
d(p* x 


"St* - 


< 1 + S 


{ !e* (r * a p* ** 5 i) + 


dy* 


r* 


t v* , 

1 + l~ r£ ) e* v* a. 




r* 


1 

L R* 




V 2 


d 

c)x* 


r* 


£S|, ^ V t N 

Sc 


sc t ) 


t y* 
1 + — — 
L R* 


dot i 
"Sx* 


+ 


+ l^ 


r* 


, . y» \( p» S . v t N. 

L R* ) \ Sc Sc J 3y* 


:]}+ W* 


(*»5) 


In the above equation Wj_ = Pe^i ^ the nondimensional mean 

rate of production of species i. The relation Vt/U e = &fc/L is used, and the 
quantity ^ where is the laminar kinematic viscosity, can be 

interpreted as the ratio of the viscous layer length scale 5 V and the defect 
layer length scale 6t : 


15 



Vt 



ik6) 


Thus, in turbulent boundary layer type flows, as a result of making the 
equations nondimensional by using the variables in Eqs. (kl) and (43), two 
small parameters, &t/ L = v t/ u e 3X1(1 £ = &v/&t appear which correspond to three 
regions of flow field: the inviscid region where the length scale is L, 

the defect -layer region where the characteristic length is St and the viscous 
layer region where the characteristic length is S v . Mellor (Ref. 10; see 
also Ref. 11 ) has shown that expansion in one parameter St/L is sufficient 
as £ = o[(&t/ L ) n ] for arbitrary n. From Eqs. (35) and ( 36 ) the nondimensional 
eddy diffusivity, v*, is: 


* * 1/2 
vt = e* ' 


l* 


(*»7) 


In nondimensional variables, Eqs. (30)-(32) and (34) are: 


x-momentum: 


d(p* u*) 
dt* 


r*“ ( 1 + _Ji Si 

' L R* 


C 1 * 6 p * u * 2 ) + I 


By* 


„ , t v*\ - ~ ~ 

r* ( 1 + p* u* v* 


t p* u* v* 


R* + y* ( 1 + 


t y* 


r* a (p* e + p* v|; 


L R* 
2 


Bp* 

Bx* 


r* 


1 + _t £ 


L R* 


6 t^ a 




5u* 

L R* 


+ 2( - 


v* 2 ** 

■= div V* 

5 *. 3 


R* + jf- y* 


A B_ 
L J By* 


^ ( 1 + J -|1 ( P *" e + P * v *) 


L R* 


°t 5v* , 

~ + 


L_ Bu* 
6 t By* 


u* 


R* + jf y* 


6 X 2 (p* t + p* v*) r 
+ 1 - 


R* + — y* 


1 1 t bv* 

\ B t y* L 

1 + L R* 


L_ 5u* 

6. By* 


U * 


R* + rr- y* 
Jj 


IT ) (p* e + B* v*) ^ 


2a , 

— u*sina + 
r* 


~ O \ ' 

+ — v*cosa) - - div V* J 


sim 


(48) 
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y- momentum: 


ft d(p* v*} 
L dt* 




{ Is (r* a 5* * v.) + w 


M , ty* 

^ U + l"R* ,X 


X p* V* 


}- 


p* U* 2 

\ 

R* + ~ y* 


L_ dP* 1 

5. 3y* „ / 6 + 


\-f a 

35E* x 


X 


r* a (p* € + p* v*) 


1 ft 5v* L_ Su* 
t / V 5 fc y* L ^oc* 6 fc dy* 

1 + L~ R* 


u* 


R* + jr- y* 


t d 


+ L 5y* 


+ T 1 IS ") (P* e + p» V*) f 2 


1 1 Tpy^ t+prv t^^-3 


^ 2 div R 


6 fc \2 (p* 6 + p* V*) 


R* + jf- y* 


2 5u* t v* 

1+ fy*^ L 6 


2 ^ 
- div V* 


L R* 


R* + y* 


r) (P* e + P* v*) 


Pp ~ ~ 2 ^ 

— ( u*sira + — v*cosa) - ^ div V* 


cosa 

(49) 


Energy: 

5 


3P 


^ a (i + r£ ) (p* H * - p*) 


+ ^ ** H*) + x 


cfy* 




r* a ( 1 + ^ g ) p* v* H* 


6 t\ 2 & 


3x* 


„a / -p* e , p* v ! 

** ( Pr~ “Pr~ ,x 


X 


1+1*22: 

L R* 


H* - 


u* 2 + ( 6 , v*/l) 2 


+ ly* 


f £*J + P * I * f I— 

\ Pr Pr^ ) by* 


H* - 


u* 2 + (6 v*/L) 2 




. , 5 * \ z a 

+ 1 L" ) 5? x 


(50) 
Contd . 
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x a r* 


( p* 6 + P* V, ) U* 


i 

L R* 


5u* _fc V* 

5k* L 5. 


R* + ^ y* 


2 

- - div V* J + 


+ (p* 6 + p* V*) j— V* 


1 t_ 5v* + L_ 5u* 

6 fc y* L 5T 6 fc W 


U* 


1 + — +- 
L R* 


R * + y* 


- - p* u* e* 


5 J „a / _^t y* 

+ l 5y* "( r * \ 1 LR* 


(p* e 4- p* v*)u* x 


1 t 3v*- L 5u * 

6 ty* L ^ ^ 


u* 


1 + _ £_ 
L R* 


R* 

Li 


u. 

+ ( p* e + p* v* )v x 


x < 2 ^-f div V.) - f ji 5» ^ 


(50) 


or in the alternative form: 
1 


d(p* h*) + 

at* 


r* a ( 1 + Jl Ft 
1 x L R* 


Is* P* **) + 


^ a (i + i£V 

1 l r* y 


x p* v* h* 


r* a (l + J2 
1 x L R* 


B t\ 2 5 


dx* 


f p* e , v fc , 
** ( — + -PJ- ,x 


5 h* 


i + l£ 

L R* 


by* 


r* a ( i + S 2*) ( £ll + flit |j* 

1 L R*/ V Pr Pr t ) 7%* 


+ ^ +_ %£^ + *^ + (^ 8(5 ** + ' P '’ t> 

1 + L R* 


1 + _t ^ 

L R* 


5u* t V* 


x 5? r 


R* + ^ y* 


t <}v* L Su* 


u* 


6 t v* L + 5 t ^ ~+ 

R* + r>* 


+ 2 I fjj ) + ^2 ( “*S ln “ + JT T* c “» : ) 2 - f 


l 5u* , s 

^~T^ + (51) 


i + 

L R* 


Contd. . 
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6 t u* sirD! + 5 fc ^ coaa: L n2 

+ l 67 — + Ss* + a F* ) 

R* + ^ y* 

x p* e* div V* ( 51 ) 


The nondimen sional equations of motion for mean turbulent boundary- 
layer type flows obtained by retaining terms of order unity and terms of 
order (S^/L) in Eqs . (44), (45) , (48)-(5l)> which include in a consistent 
way longitudinal curvature and normal pressure gradient, are: 



Continuity: 




(1 + — ^ 

1 L R* ) 


8 / „a ~ N , 8 

(r* p* u*) + 


8x* 


cp 




p* V* 


r° 

(52) 


Species Continuity: 


B(p* 02) 
Bt* - 


r*~ (i + 3 SI 

1 L R* 


|e* (r * & ~ p * * S 1> + [ 


r» a (l + Jl! ) X 

1 L R* ; 


X p * V * oc 


1 8 

a dy* 


r* 


r* 


1 L R* ) \ Sc Sc t J ciy* 


, T \* 
H- W. 
1 


(53) 


x- momentum: 

8u* , u* 


3* 


bt* 


8u* , ~ 8u* ( t U* V* 


1 + 

L R* 


+ v* 


L R* 


1 + _t 22 

L R* 


_1_ B_ 
a by* 


r* 


r* 


5*vt> (S-iSg 


8u* t U* 


6 fc (p* e + P* v|) s ~ # 


R* 


By* 


(5^) 
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y- momentum: 


-* / 3v* „ 

^ p l sf* + u 


3v* ~ 3v* \ p* U * 2 _ L 39* °t 1 3 

Sc* "Sy* ) ~ R* 6 fc Sf* L r *a ^ 


Si / — A “ ■X" \ Su* 

r* (p* € + p* v fc ) ^ 


_t _i_ 3_ 

+ L 3y* 


r* a (‘p* £ + p* v*) (z^£ 


2 & 
■5 dlv V* 


(55) 


Energy: 


3_ 

3t* 


+ irS ) (p* s * - p*) 


+ ^ (r* a p* u* H*) + ^ 


dy* 


r* a ( L + _t sr ]x 

1 1 l r* ’ 


x p* -v* H* 


3 

3y* 


r* 


c j* V e P* 3 f ~ u* 2 

1 + rpjlir + ^- — 


3y* 


,a/, . 6 t y* > a , ^ ( 3(u* 2 /g) 5 t u* 2 ^ 

R* 




(56) 


or 


P* 


dh* , u* 6h* -v Sh* _ 1 B 

sr + — 67^55 + ^5s;=^:^ 

1 + L R* 


^f 1+ ^^V^ + 

' L R* J V Pr 


p* v £ \ 3h* 


Pr t JW 


3p* u* 3p* ~ 3p* ,-a - v 

5t* 5 ^ 3x* + "Sy* + ^ e + p* v*) 

1 + L~ R* 


f 3u* 'N 2 _^t 3u* u* 

VSPy " l 3y* R* ■ 


(57) 


4. TWO- EQUATION TURBULENCE MODEL FOR CURVED BOUNDARY LAYER TYPE FLOWS 


The turbulence quantity defined in Eq. (26) and satisfying the model 
equation, Eq. (36), is regarded as the total kinetic energy of turbulence, with 
u" and v" denoting the fluctuating velocity components in the streamwise and 
in the normal to the shear plane directions, respectively (in this section we 
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consider thin shear layers parallel to y- constant lines) . However , Wilcox 
and Chambers (Ref . 2) argue that this definition of e is physically realistic 
only if the kinetic energy of turbulence is equip ar t it i one d, i.e., if the 
turbulence is isotropic. Stating that turbulence in boundary layer type 
flows is anisotropic, they postulate that, more appropriately, e is propor- 
tional to v n 2. The proportionality coefficient is then established by 
requiring that e be numerically equal to the kinetic energy of turbulence in 
the law- of -the- wall region (Ref. 5) of a flat-plate boundary layer, i.e., 
e = 9v n '^/4 . Furthermore, by proper consideration of the physics of turbulent 
flows with streamline curvature based on the classical stability arguments for 
flow over a curved wall (see Ref. 12), Wilcox and Chambers show (Ref. 2) that 
the stability of such a flow depends mainly on the behaviour of vertically 
moving fluid particles and not on any attendant nonvertical fluctuations. 

Hence they conclude that in a turbulent boundary layer, velocity fluctuations 
normal to the wall, v" 2 , play an important role in curved streamline turbulent 
flows . 

The model equation for the "mixing" energy, e = 9v"^/4, may be 
obtained from the exact normal Reynolds stress equation (A33) derived in 
Appendix A, rewritten in the following form: 


b f - pv f y j 

St 


r a (l+Ky) 




a - ~ pv’V r 
r p u 


Sy 


r (l + Ky)p v 


■nrrrr 


~ P v 7 
P 


]} 


2k 

1+Ky 


u "pu TT v Tr 


= - 2pu rr v 



1 dv _ KU 
1+Ky c*x 1+sy 


pv"v" 2 


^ ■ 


Production 


2r *y _ 2 + 

l+*y yy dy 

Dissipation 


Q — { [r a (2r v" - u"pv"v” ) ] + 

r a (l + .y) l^ L ^ *3T 

+ ^ [r a (l + Ky)(2T yy v" - v"pv'V‘ - 2p V') ]J- + 


— Laminar, Turbulent and Pressure Diffusion 


(Contd. . . ) 
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+ + 7717 v"(p' - Txx ) + |S v'V - T^Jcosa + P« 


l+hty 


l+pcy 


W 


"5y" 


Turbulent Work 


-^l 


( 58 ) 


Compressibility 


Multiplying Eq. (58) by 9/4, modelling the groups of terms for the 
corresponding physical quantities as for Eq. (A3*0 and taking into account 
Eq. (28), we obtain the f ollowing model equation for the mass-averaged 
"mixing” energy e = 9v”v”/4p: 


S(p e) , 1 f £(r a p u e) , 5 r a /-, . \- ~ -,1 , 

^ 7777 I 5y [r (1 + v e] j + 


+ 2 * u " ( 1 5v 9u KU 

2 1+icy pv t ( 1+ny c5x )§y 1+ny 


x? 1 r ^ r r6 Ct* + °* p v t ) Se 

(a*ps - p* e)e - |* e p div V + Q 1 ^ 1 1 de 


r (l+ny) 


1+icy 


a 


+ |f 


r a (l + Ky)(|j, + o* p v t ) || 


(59) 


Written in nondimensional variables of Eq. (4l) , Eq. (50) becomes: 




a f n t # 
r* 1 + — 1 

1 X L R* 


S(r* a p* u* e*) 5 

— s# — + a* 


+ rp)^"* e * 


+ § 


u* - v * 

T P* V 4- 


R* + y* 


1 °t Sv* ^ L 5u* U* 


% y* L ^ 5 t ^ ’ R * 
1 + L R* 


= (cn* p* S* - (3* e*)e* - l e* p* div V* + 


x 


r* a f 1 + — — 

1 L R* 


(60) 

Contd. 
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f( & t \ 2 h 

- r*( p* e + cr» p* v*) ^ - 

x lr ) ssf 

555 J 


L , R* 


+ 


b_ 

By* 



l 


6, „ 

t y* 

L R* 


(p* e + a* p* v*) 



(60) 


where 


s* 



i + 2 l 2* 

- L ‘ ' -r- T-iJV 


L R* 


\ Bk* 


L R* 




1 + ^ 


t y* 


L R* 


^t &y* 
L chc* 


L_ du* u* 
5 fc 3y* " R* 


-P 


(61) 


get: 


Retaining terms of order not higher than 6^/L in Eqs . (60)-(6l) we 


d(p* e») + 

5t 


r* a ( 1 + 


_1 f B(r* p*u* e*) B 

6 t \ SE* dy* 

l R*y 


** a ( 1+ rl e * 


Ju*- „ L 6u* 
2 R* P* v t 6 fc dy* 


a * 5 *%^-Si- p * et 


e* - 


- |* e* p* div V* + 


1 


- a (^S) 


dy* 




x (p* e + cr* p* v*) |p | 


(62) 


Except for the compressibility term, which is retained here in its 
general fom, Eq. (62) reduces to the corresponding equation in Ref. 2 for 
a = k = 0. 


Writing the dissipation rate equation, Eq. (37) * in nondimen si onal 
variables, Eq. (4l), we have: 
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d(p* e* ) . 1 f 5 , ,a ^2 N , 

Va 7 , p * u *^ ) * 


r* 1 + 


L R* 


+ h* P* C 1 + IT 8 ) " P * ^ 1 } = a “ P * ** “ 


P + 2cr 


i 6 t a.e* 'n 2 , /" a^* \ 2 1 1 ^3 , 

—57^-^; + ^J Jj 6 * 

1 + L R* 




6 tf a r r* (p* € + a p* v t ) 5e *2 


L R* 


** l 1 + rS ) ( _p *^ + CTp * v V^r 


. _^2 ,. .5* e* t div V* f ap* ap* sp* dp* A , 

- I 2 p. s* 4LVV. - 5 3 ^j ^VlS-cF-^S*; (63) 

1 + L R* 


Retaining terms of order not higher than 5^/L in Eqs . ( 63 ) we get: 


5(p* e* 




a(r* a p* u* e* 2 ' 


( 1 + r S ) "p* ^ e * 2 


a P * fc- 


>L au* u* 


6, ay* r* 


p + 20- 


e * 2 + 


W i + II 2 1 

V L R* 


b C t *“ ( 1 + r § ) 6 + 1 _c * V P l£- 


1 + L R* 


2 4 



Equations (62) and (64) represent the model equations for the 
quantities e* and e* for boundary- layer type flows (&t/L « l) needed to 
determine the turbulent (eddy) diffusivity v*. The values of the closure 
constants j a, a*, 0, 0*, cr and a* and |* 5 are the same as in Eq. (28) . The 
magnitudes of the coefficients £2 and S 3 and hence the importance of the 
terms which contain these coefficients may be determined by comparing 
theoretical predictions with reliable experimental results. 
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APPENDIX A 


DERIVATION OF TWO-EQUATION TURBULENCE MODEL 
FOR COMPRESSIBLE PLANE OR AXISYMMETRIC TURBULENT FLOWS 
IN ORTHOGONAL CURVILINEAR BODY-ORIENTED COORDINATES 


The mass-averaged turbulent kinetic (mixing) energy e, and the 
mass-averaged turbulent dissipation rate €, are needed to define the eddy 
diffusivity given by Eq. (35) onp.32. The following detailed derivation 
of the model partial differential transport equations in the coordinate 
system considered which e and e satisfy follows the Wilcox and Traci approach 
(Ref. 5) of turbulence modelling. 


1 . The Turbulent Dissipation Rate Equation 

^ Tfte dissipation rate quantity is chosen as S = pCi[ Q^/p where 
£l± = % + %, ^i = pu>i and the vector a? = rot v* is the vorticity. In vector 
form the momentum equation may be written as (see, for example, Ref. 7, Ch. 
1 ): 




+ ^ 


2 


+ 


-> 

0) 



(- p ? +% 


(Al) 


where I and t are the unit and viscous stress tensors respectively. Taking 
the rot of Eq. (Al) we get, using the vector relationship rot (a? x v^ = 

'=* ^)S? - (u • , 




St ■ 


(v 


V)i 




V^v + ■ v) = rot 


^ • (■ 


I + t) 


(A2) 


From the continuity equation (l) we have ^ • v = - (l/p)(dp/dt)and Su^St + 
(t • = d£?/dt. Hence Eq. (A2) becomes 


P 



(fM • - 


dp 
W dt 


p rot 


4 =$ & 

‘ (- p I + t) 
P 


or 



^)v + 2f?(V^ • v) = p rot 



(- p I + t) 


(A3) 


The right-hand side of Eq. (A3) may be written as 




(- p I + t) = _ gradp 
P P 




T 


P 
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The components of the vector V • r are* in the orthogonal curvi- 
linear body-oriented, coordinate system* considered: 


x- component: 


0 x a. 


r (1+Ky) 
y- component: 


7— ' { s < r \*> + |f + V 1 } + ts? - ? T W S1 ”“ 


e.. = 


y a 


r (l+Ky) 


\ (r a T ) + [r a (l + Ky)r ]T - - — t cosQ ! 

dx v xy' oy L v yy J J 1+xy r cpcp 


and 

.^=n a + n^+ 0 • ~y 

P P P 1 x y p 7 

where C? 5 p* and 7* are the unit coordinate vectors in the x* y and cp directions 
respectively. Hence 


a 


7 


rot 


gradp ^ = ^2 S ra(i P x gradp = — 
P O 2 


1 dp 

h 2 ciy 

1 5p 1 5p 


h 1 "Sc 


h 2 dy 


^ 0 


(l+Ky)p‘ 


/ dp dp dp dp \ 

2 ) 7 


and similarly 


rot 


^ • f 


p J ' l+Ky { Sx 

Equation (A3) may thus be written as 
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- L 7 


"SF 


p ^ - p(j? • ^)v - 2pi?(^ • v) + ^ + { 


an 

y 
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a[(l+xy)n x ] ^ 
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(a4) 
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For the only component ^3 = & of the vector S? we get from Eq. (a 4) , 
in the coordinate system considered, 

cto , ou ^ _ cto « n r1 . Tr . 1 f Sp 5p Sp 5p \ . 

p 5t + I+^F ^ + pv " " 2pfi dlv 1^ ^ ^ S J + 

Decomposing ft = ft + ft" , using Eqs . ( 19) and. (20), and multiplying 
Eq. (A5) by r a (l+Ky)ft” and taking the average, we obtain, term by term: 


P&” g£ [r a (l + Ky)(n h- ft")] - ^ g£ r a (l + Ky)p ~ 


a /n ^ Q'to" ^p 

r (1 + Ky) — ^ 


(A6) 


r a pfi"(u + u”) J | 


2 Sc 


p u + u"'pn"ft"' 

P 


FF b 

2 Sc 


(r a pu) + 


r a 


bn 


(A7) 


r a (i + Ky)pn"(v + v") (n + a") = \ ^ r a (i + k: 


y) 


- ~ p n fl _ + v 't pn <t n .> 
p 


} 


*Mj- [r a (l + Ky)pv] + r a (l + ny) pv"£i" |p 


(A8) 


- 2r a (l + K y)pF(fi + o")div(V + V*) = - 2r a (l + K y)p div V - 

P 


- 2r a (l + K y)[H pfi" div V*' + pn"F div V*'] 


(A9) 
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5(p l± p { ) . 5(p + p t ) _ S(p + p' ) . 5(p + p' ) 
dx dy dy 5c 


r a a" 


= r 


771T / 3g. 3 p 3p dp \ /sit 5p' . cp ^11 3 p' 3p 




(A10) 


and the last term in Eq. (A5) nay be written as 

an. 


0 s jV - 5 CCi+-y>Ii t l} 


(All) 


Adding Eqs . (a 6)-(A11), taking into account the continuity equation, 
Eq. (ll), and regrouping the terms we ge t the following equation for the 
turbulent dissipation rate quantity e 2 = pft ,T ^”/p; 


h e2) + 1 


r a (l+Ky) 


^ (r a p u e 2 ) + [r a (l + Ky)p v e 2 ]j- = 


pu dn . n^-vTt 

r^T 5E pv n dy 


Production 


(r a ^W) + ^ [r a (l + Ky)v”pn"n"]j- 


r a (l+Ky) 


Turbulent Diffusion 


l+Ky \ cix §y j 


— Laminar Diffusion and Dissipation — 
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o ~ ~ 

- 4p € div V - pft" div V" 


4pfi"fi" div V" + 


1+Ky 


dp dp _ dp dp 
cfc ^ ” c5y "Sx 


+ 


+ n 


it 


5p* 5p* _ Sp T dp 1 \ 
dx w~ ’ dy dir J 


+ ft" 


dp 1 dp 
ciy~ Sc 


+ n" 


dp 1 dp 

ST "5y 


n" 


dp T dp 

ST dx 


fl" 


dp 1 dp 

ST dy 


Compressibility 

(A12) 

The physical meaning of each group of terns is indicated in Eq. (A12) . Each 
group of terms should be "modelled" to express the physical quantities actually 
present (rates of turbulence production, dissipation and diffusion of the 
quantity e 2 as well as the effect of compressibility on the rate of change of 
e 2 ) in functional forms containing a closed set of dependent variables. These 
functional forms should be dimensionally consistent and may contain empirical 
constants that are expected to be insensitive to the character of individual 
flow fields. The following groups of terms are modelled after Wilcox and 
Traci (Ref. 5): 

Production: 


Using the eddy diffusivity approximation, 


_ u "^ n _ pV t cto 

pu ~ 1+^ SE 5 


pv ■ pv t IP 


this group of terms may be written as 


2pv, 


7 i bn 'f 

. ( f 1 

\ i+ny dx J 

+ Vay J 


(A13) 


and is modelled by Wilcox and Traci by the dimensionally consistent expression 


a p S e 


(Al4) 


where 


S = 2(e e + e e + 2e e T 

Y xx xx yy yy xy yx' 


(A15) 
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_ - Sv 

yy ” ciy 5 
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xx 


1 

l+Ky 


5u 

3x 


+ KV 


_ 1 f _1 3v + 5u _ KU \ 

e xy e yx 2 ^ 1+sy ’§x dy 1+Ky ) 


(A16) 


and a is a closure coefficient. 

Dissipation : 

The group of terms describing the dissipation is modelled as 
follows: from dimensional consideration the rate of dissipation may be 

taken as e/p[l/sec] and the dissipation rate of the quantity p e 2 in Eq. 
(A13) may be taken as 


- P P e 


2 



(A17) 


where p is the proportionality coefficient. 


Laminar and Turbulent Diffusion: 


The form of this group of t erms and, the eddy diffusivity approxi- 
mation for the quantities u ,, pft ,, G n and v n pft ,r ft lT suggest the following modelling 
of these terms: 


r a (l+Ky) 


5x 


a / 

r r (n 


+ cr pv t ) 


1+Ky 


Sx 


+ b 


2 


r (1 + ny)(p. + cr pv t ) 


where c is a closure coefficient. 
Compressibility : 

The first group of terms 


2 


- 4(p e div V + ft pft" div V M + pft"ft n div V") 


- 2 

is missing from the corresponding equation for p e written in general 
three-dimensional' Cartesian coordinates in Ref. 5- This term is modelled 
here simply as 


- Sc 


- 2 

P € 


73 

div V 


(A19) 
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where ^2 is a proportionality (closure coefficient) . It is worth noting 
that this term is similar in form to the compressibility term for the 
turbulent kinetic energy in Ref. 5 or in Eq. (a4o) . 

The second group of terms in Eq. (A12) labelled "compressibility” 
is neglected in Ref. 5. However, this term may be retained and modelled 
in the following way: from mass- averaging relationships we have 


fi" = _ (A20) 

P 

Consider now the idealized situation of one -dimensional high-speed flow for 
which the total temperature is approximately constant, the fluid is a perfect 
gas and pressure and total temperature fluctuations can be neglected. Then 
from the energy equation 


7-1 2 2 
-hr- u + a 


7+1 

2 


a* 


where a is the speed of sound and a* the critical speed of sound. With 
u = u + u* and a = a + a f we get 


(u 2 + 2 uu* + u’ 2 ) + a 2 + 2 aa* + a * 2 = ^ — a * 2 
(constant total temperature) 

-2 -2 2 

Since ( 7 -l)u /2 + a = ( 7 +l)a* /2 we have 

( 7 -l)uu’ + U' 2 + 2aa' + a ' 2 = 0 (A21) 

2 - - 

On the other hand pa = 7 P; with p = p + p 1 and a = a + a* (pressure fluc- 
tuations are neglected) we get 


— (a 2 + 2aa T ) = - 2aa’ + a' 2 
P 


or 


P 1 


2a 1 + 3a T 


-2 

a 


(A22) 


Equation (A2l) gives 


a a a a 

Substituting the value of -2a* /a from Eq. (A22) into the above equation, we get 

uu f _ p l 


(7-1) - 3 1 

a p 
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or approximately 


P T _ 
P 


( 7 - 1 ) 


(A23) 


For two-dimensional flows we can generalize this equation to 


— ^ (uu t + vv’ ) 
p a 

Using Eg. (A24), Eg. (A20) may be written as 


(A24) 


(u u’n 1 + v v f 0,' ) 
a 


or using again the eddy diffusivity approximation 


- u T ft' 


v t be 


1+Ky J 


- v’n 1 


be 
v t by 


we obtain 


. r v t f u be , - be \ 
~ -2 \ 1+Ky 


(A25) 


However, from Eg. (A12) for e , in the approximation considered, we can write, 

u be ~ be .7 

l+Ky Sx dy 

Taking into account Eg. (35) , Eg, (A25) becomes 


- In ^2 e div V = 
^ a 


i 3 € dlv V 
a e 


lo div V 
J a 


(A26) 


where £3 is a closure coefficient, 
terms in Eg. (A13) is modelled as 


Hence the second group of compressibility 


which is dimensionally sound and similar in form to Eq. (A19) . From Eqs . 
(Al4) 5 (A17)-(A19) and (A20) , the resulting model equation (A12) for the 
rate of change of the quantity e 2 , is then 


Mp £ 2 ) + i 

dt „a 


by 


r (l+icy) 
= a p s e 

2 i 


^ (r a p u e 2 ) + [r a (l+Ky)p v e 2 ]j- = 


-{ 


£ + 2o- 


1 


e 3 + 1 


r a (l+Ky) ^ ^ 


1+Ky Sx 

r r& (p + P v t ) ^ € 2 

1+xy dx 


by 


be 


r (1 + xy)(n + o- pv t ) ^r 


„ - _2 „ p e div V f dp $p cip 8 p 'N 

- i 2 p € iiv V - £3 “ 1+^T ^ 3y - ^ ^ ) 


(A28) 


The dissipation-rate equation (A28) contains a term proportional to (l/l+rcy) 
x(cl#/5x ) 2 + (<3$/8y ) 2 referred to as the gradient- dissipation term in Ref. 5 • 
Its introduction was motivated by comparison of the Wilcox-Traci model with 
other two-equation turbulence models, particularly those developed by Ng and 
Spalding (Ref. 8 ) and by Jones and Launder (Ref. 9 ). It enables one to 
predict more realistic values for the turbulence length scale especially 
near solid boundaries. The values of the closure coefficients a, f3 and cr 
are given by Wilcox and Traci (Ref. 5) as: 


2 x 



1 - ±° 
x 11 


ejqo 



Re, 


Be 1 / 2 ,. pe 1 / 2 

H 3 1 - — 

(A29) 


The magnitudes of the coefficients £2 £3 ^y be assessed by comparing 

theoretical predictions with reliable experimental results. 


2 . The Turbulent Kinetic Energy Equation 

The x-momentum equation (13) ma y be written as 


r a (l + ^y)p ^ + r a pu + r a (l + «y) pv ^ + r^Kpuv = - r 


Su _ L „a ( 
8 x 


8 u T _a 
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a 8 p 
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dx dy L ' ' xy J xy cpcp dx 


(A30) 
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Let 


u=u + u", v = v+v", p = p + p’ , P=P+P ! 


and 


t, . = r. . + t! . 
ij ij 13 


(A31) 


Multiply Eg.. (A30) by u” and take the average j after some calculation we get 
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or rearranging the terms, 
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written as 


Using Eq. (A3l) and multiplying the y-momentum equation, Eq. (l4). 
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Combining Eqs. (A32) and (A33) we obtain for the mass-averaged turbulent kinetic 
energy, Eq. (26) , the following equation: 

^ + (pa " p s e) + If [ra(1 + Ky)p 7 e] } = 
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Each group of terms expressing the physical quantity present is "modelled 11 in 
the following way (Ref. 5): 

Production 

Equation (28) suggests a modelling of this quantity in the form 

a* p s e (A35) 

where 


S = /2(e e +~e e + 2e e J 
I K xx xx yy yy xy yx' 

and a* is a closure coefficient. The expressions for e . 
in Eq. (Al6) . ^ 

Dissipation 

The rate of dissipation may he taken as e/p [l/sec], hence the 
dissipation of the quantity pe may he modelled as 

- P*pe | = - p*ee (A37) 


(A36) 

e . e are given 
yy 5 xy e 


where is a proportionality (closure) coefficient. 

Laminar and Turbulent Diffusion 

The form of this group of terns suggests the following modelling: 
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with o* as a closure coefficient. 
Compressibility 

Using the relationships 
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and Eq. (A24) ? we can write this group of terms as 
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From the continuity equation div(p V) = 0 we have 


Hence the compressibility term is modelled as: 


i - ^ I - - 5*«5divv 


(A39) 


where |* is a closure coefficient. 

The resulting model equation for the mass- averaged turbulent kinetic 
energy e is then 

p 5 e ) + b + Ky)_p ^ e] } = 

(a* p s - p* e)e + 
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1 Sx 
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}- 


I* e p div V 


(a4o) 


The values of closure coefficients (see Ref. 5) are: 



1 _ ii e3q> ( - 




9 

100 ’ 




(a4i) 


The turbulent Reynolds number Ret is given in Eq. (A29) . The turbulent work 
terms in Eq. (A34) are neglected. 

Equations (a4o) and (A28) represent a two-equation turbulence model 
for compressible plane or axisymmetric Navier-Stokes flows in the coordinate 
system considered. 
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